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Abstract 

The isotropy in the ultra high energy cosmic ray (UHECR) flux observed by Yakutsk 
and AGASA experiments, is a very strong constraint to production and propagation 
models alike. Most of the scenarios proposed in the literature should produce a sizable 
anisotropy as either extragalactic luminous or dark matter is normally associated with 
the invoked particle sources. We explore the possibility that the magnet ic fields in 



fossil cocoons of former radio galaxies - so called radio ghosts (Enfilin, 199E)- are able 
to scatter UHECR in the intergalactic medium giving rise to the observed isotropy. 
We show, through numerical simulations, under which conditions this process can be 
operative and the magnitude of the effect. We further demonstrate, that if radio ghosts 
mix with the ambient medium, they might be able to produce the observed magnetic 
fields in clusters of galaxies. In the case of mixing, the UHECR isotropization would be 
even stronger than in our conservative estimates. 
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1 Introduction 



1.1 The Problem of UHECR Isotropy 

The upper end of the cosmic ray spectrum, at total energies above the (as ye t unobserved) 
Greisen-Zatsep in-Kutzmin (GZK) cut-off (~ 4 x 10 19 eV - ( Greisen, 1966 ), ( Zatsepin and 
Kuzmin, 1966| )) represents a challenge for particle-physicists and astrophysicists alike. The 
nature of the sources of these ultra-high energy cosmic rays and their distance scale are still 
unknown. Only our own galactic disk can be ruled out at present as a major source site, as 
a compatible anisotropy has not been observed by any of the experiments sensitive to the 
UHECR energy range ( JTakeda, 199S| ; [Bird et al., 1999| ; pVledina Tanco and Watson, 1999[ ). 

Although photons, neutrinos, or some unknown particle cannot be disregarded, the muon 
to electron ratio measured for extensive air showers (Halzcn et al., 1995; Hayashida, 1996) 
points to hadrons as the primaries hitting the upper atmosphere and triggering the cascades. 
Neutrons with relativistic factor 7 ~ 10 11 decay into protons after a path of ~ 1 Mpc. Heavy 
nuclei, on the other hand, may lose as much as ~ 2 — 4 nucleons per Mpc (Puget et al 



1976; Cronin, 1992) due to photodisintegration in interactions with the cosmic microwave 
and infrared backgrounds. Therefore, although new estimates of the infrared background 
(Malkan and Stecker, 1998) seem to indicate that the previous photodisintegration rates may 
be overestimated by a factor of 10 (Btecker, 1998), UHECR, unless galactic, are very likely 
light nuclei, probably mainly protons. We will assume the latter in the remaining of this 
work. 

Depending on the large scale co nfiguration of the intergalactic magnetic field (IGMF), field 
values as low as Bjgmf ~ 1CU 9 G ( Kronbcrg, 1994a ) may be expected. Consequently, proton 
gyroradius at E ~ 10 20 eV can be of the order of 10 2 Mpc. This represents a potentially 
interesting perspective for astronomy since UHECR should, under favorable conditions, point 
to their sources opening a new astronomical window. This would allow not only the direct 
study of their sources, but also the measurement of the intervening intergalactic and galactic 
(particularly halo) magnetic fields. 

Relevant to this potential is the nature of the production mechanism involved. So far, both 
bottom-up and top-down UHECR production mechanisms have been devised and constitute 
perfectly viable scenarios. The first set of models invoke some kind of acceleration agent 
that is always spatially associated with baryonic matter, closely traced by luminous matter. 
Top-down mechanisms, on the other hand, resort to the decay of relics from GUT phase 
transitions in the early Universe (Bhattacharjee and Sigl, 2000). Under the assumption 
that relics aggregate like dark matter (DM) and that their decay products are susceptible of 
inelastic interactions with the CMB, then any observed anisotropy should trace the nearby 
distribution of DM. Therefore, the UHECR flux should have two components, one originated 
inside the galactic halo and another extragalactic still roughly traced by luminous matter. 
(Dubovsky and Tinyakov, 1995) showed that, under general conditions, the halo component 
would dominate the extragalactic flux by at least two orders of magnitude. This is only true, 
however, in the unrealistic case of dark matter uniformly distributed in intergalactic space. 
Nevertheless, dark matter aggregates strongly and tends to be overabundant, by factors of 
~ 10 2 , in the center of galaxy clusters wh en compared to its abun dance in the halos of isolated 
galaxies. It can therefore be shown that ( Medina Tanco, 2000a ), in a sample of 47 events (as 
AGASA's above 4 x 10 19 eV), and assuming Virgo as the only source of extragalactic events, 
3-7 events should originate in Virgo and arrive inside a solid angle of approximately the size 
of the cluster. This could give rise to a slight anisotropy that correlates with the SGP when 
combined with the almost isotropic flux originated in a large galactic halo. 

Photo-pion production interactions with the cosmic microwave backgrou nd (CMB) re - 
strict the source region to nominally D < 100 Mpc (Waxman et al., 1997; Hillas, 1995). 



Therefore, the signatures of the characteristic inhomogeneity of the local universe should be 
clearly observable in the angular arrival distribution of UHECR regardless of the production 
mechanis m. In spite of the low statistics availa ble at present, the latter does not seem to be 
the case ( fTakeda, 1999| ; |Medina Tanco, 2000b[ ). 
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If UHECR are charged particles originated inside the large scale structures in the nearby 
Universe, then any isotropization must be due to intervening magnetic fields. This could 
happen cither inside the galactic halo, like in the presence of a magnetized galactic wind 
( Aim ct al., 2000|), by the in teraction with highly structured IGMF inside walls and fila- 



ments ( Medina Tanco, 1998b| ) or by the scattering off magnetic irregularities permeating the 



intergalactic medium. 



In the present work we analyze the latter possibility considering radio ghosts (Enfilin 



1999) as scattering centers of UHECR protons. 



1.2 The Existence of Radio Ghosts 

Active galaxies eject large amounts of radio emitting plasma - short: radio plasma — into their 
environment. There it forms the typical cocoons of radio galaxies. The radio emission results 
from synchrotron emission of a population of relativistic electrons in the radio plasma's 
magnetic fields. It is possible to estimate the minimal energy density of the electron and 
magnetic components of the radio plasma required in oder to produce the observed emissivity. 
This minimum is given by rough energy equipartition between electron population and the 
magnetic fields. The resulting minimal pressure is typically of the order of the environmental 
thermal pressure, indicating that relatively strong magnetic fields are present. 

After a cosmological short time of 10 7 — 10 s years the radio luminosity of the cocoon 
decays strongly due to radiation and expansion energy losses of the electron population, or 
since the central engine of the radio galaxy stopped its activity. Although undetectable by 
our instruments, the radio plasma is still present in the IGM inside a fossil radio cocoon, 
a so called radio ghost. The subsequent evolution of the radio plasma is unclear, since 
observationally poorly constrained. It can be expected that the strong magnetic fields of the 
ghosts allow it to resist erosion by subsonic turbulence. The fossil radio cocoon should be 
kinematically decoupled from the ballistic motion of the parent galaxy and follow mostly the 
flow pattern of the embedding material. From this one would expect the ghosts to have a 
cosmological distribution comparable to that of the galaxies. And one would further expect 
that the oldest ghosts are swept into clusters of galaxies by the flow of structure formation. 
But buoyancy of the probably very light radio plasma can produce some relative motion 
between ghosts and the IGM gas. That buoyancy is operating on radio plasma is impressively 
demonstrated by a recent simultion of a buoyant rising radio blob in a cluster center, which 



reproduces radio and X-ray observations of M87 strikingly well ( Churazov et al. , 200C ) . 

This might allow the ghosts in clusters or filaments of galaxies to ascend to larger radii, 
until they get stopped by freely infalling matter at the accretion shock. It depends crucially 
on the topology of this accretion shock surface if the ghosts are able to escape from gravita- 
tionally bound structures as clusters and filaments of galaxies, or not. It is therefore difficult 
to predict the spatial distribution of ghosts. In this work we consider cases where the ghosts 
are distributed as the galaxies, where they are more, and where they are less clustered. 

There are several observations supporting their existence and future experiments might 
be able to see further signatures of ghosts. 

• The progenitor of radio ghosts, radio cocoons of active radio galaxies, are observed and 
well studied. Since the sharp boundaries of radio cocoons indicate that these objects 
are able to resists their very turbulent birth, it is reasonable to assume that they also 
survive the much weaker ambient turbulence. 

• A radio ghost can be revived to radio emission by the compression in a shock wave. 
This should lead to regions of diffuse synchrotron emission within clusters of galaxies, 
where accretion or merger shock waves are expected. Such objects are indeed observed, 
and are named cluster radio relics. Models for such processes can be found in Enfilin 



et al. (1998a) and Enfilin & Krishna (2000). 



• Another class of diffuse radio emitter in galaxy cluster might also be connected to 
ghosts: cluster radio halos. There is the speculative possibility that the short-lived 



3 



radiating electrons are replenished by the decay of charged pions produced in hadronic 
collisions of long-lived cosmic ray protons, which were released from their confinement 



in radio ghosts (Enfilin, 199S) 



• Although at low energies, the electron population within ghosts should remain mostly 
relativistic. Comptonization of the CMB photons by this electron population should 
produce CMB distortions, which have a spectral signature comparable to th at of a grey 



absorber. Th is signature is weak with an optical depth of the order of 10 7 ( Enfilin and 



Kaiser, 2000 ), but principally observable. The energetically up-scattered CMB photons 



might also be detectable. 

For this work, to estimate the deflection of UHECR by ghosts, most of their above discussed 
properties are not of any importance. Only the magnetic fields are required, and their 
existence is nearly guaranteed by the existence of radio galaxies, and the extremely low 
magnetic diffusivity in extragalactic plasmas. Therefore, in the context of this work, the 
term 'radio ghost' can be simply translated into 'remnant magnetic fields of a former radio 
cocoon'. Even if turbulence is able to disrupt ghosts, the shear motions should increase the 
magnetic energy content of the fields, leading to an enhanced efficiency in UHECR deflection. 
We summaries, that the existence of regions in the IGM filled with remnant magnetic fields 
of former radio galaxies - as it is assumed in this work - is difficult to be doubted. 



2 Properties of Radio Ghosts 

2.1 The Distribution Function of Radio Ghosts 

The properties of a radio ghost depend on three factors: the state of the radio plasma, as 
released by a radio galaxy (RG), the turbulent history of the environment, and the present 
environment of the ghost. We are mainly interested in ghosts located in filaments of galaxies, 
where turbulence is less strong compared to galaxy clusters. In our simplistic approach, we 
neglect therefore any environmentally induced changes in the ghost properties. We assume 
that a ghost was originally produced with a total energy P g h, and this energy did not change 
significantly over cosmological times. The ghost consists of magnetic fields and relativistic 
particles (electrons and protons), which should be roughly in energy equipartition. The 
ghost's internal energy density 

£gh = £b + £cr ~ 2 Eb (1) 

is related to it's internal pressure via the relativistic equation of state, which the CR popu- 
lation and the tangled magnetic fields do follow: 

1 2 

fgh = g^gh ~ • (2) 

Since the ghost is in pressure equilibrium with its environment P g h = P, the energy released 
by the RG can be written as the sum of the internal energy of the ghost and the volume 
work required to produce a cavity in the IGM for the ghost: 

E gh = e gh V gh + PV = 4PV gh . (3) 

The relevant properties of a ghost, as volume Vgh and magnetic field strength P g h, are 
therefore determined by P and E^: 

V gh (£ gh) P) = || (4) 

P g h(P) = V12ttP (5) 

Thus we need a model for the IGM pressure. The enthalpy of the gas is increased during 
structure formation by compression of the flow, which is driven by the evolving gravitational 
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potential]^] . This picture neglects non-gravitational heating and cooling processes. Cooling 
is only important in the center of dense clusters of galaxies. Whenever the gas is subsonic, 
the pressure is related to the gas density via 

P(n gas ) = Po (%as/«gas,o) 7 , (8) 

where 7 = 5/3 is the adiabatic index of the gas, otherwise the pressure is lower. We neglect 
deviations from this equations assuming that only a small fraction of gas has a sonic or higher 
velocity. This is justified, since supersonic gas becomes soon shocked when it hits the next 
accretion shock around a cluster or filament of galaxies. 

Unfortunately, the distribution of gas in the local Universe is not observed with current 
instruments. But the number density n ga i of the galaxy distribution can be used to give a 
rough idea of the gas and pressure distribution. We assume that the baryon fraction of the 
matter density is approximately the same in all structures (f? gas /f?matter = const), and that 
the bias between optical galaxies and matter can be approximated by 

^-gal ^-gal,o(^'gas/^'gas,o) g • (9) 

The galaxy bias can be found to be 6 ga i = 1.6 by lin earizing this ex pression and comparison 
with the observed linear bias for optical galaxies ( Peacock, 1999| ). Thus we get a rough 



proportionality between galaxy density and pressure: 

P(n gal ) = P (n ga i/n gaLo )^ 6 -' ~ n^ 42 . (10) 

This relation can be normalized by using the observed properties of clusters of galaxies. A 
typical density and temperature of a cluster is n e = 10 -3 cm -3 and kT = 5 keV, leading 
to P a = -Pciustcr = lOkeVcm -3 . We therefore get a ghost magnetic field strength of 25 /iG 
inside of galaxy clusters and 0.54 /iG in galaxy filaments (assuming n e = 10 -5 cm ). 
The distribution of ghost energies is assumed to follow a Schechter distribution 

n gh (E gh ) dE gh = N gh (§^) exp (-^)^, (11) 

for E g h > Ejnin, ngh(-Egh) = otherwise. Both the normalization factor iV g h and the upper 
cutoff E* are allowed to depend on the position, parameterized by the galaxy density (iV g h = 
Ngh(n g ai) and E* = E*(n ga i)). 

It is convenient to introduce the ghost-volume-filling factor 

POO 

$ gh = dE gh V gh (E gh ,P)n gh (E gh ,n ga i), (12) 
and the average energy density in ghosts 

4 

e g h = g £ g h$gh. (13) 

The factor | appears here, since our definition of the energy of ghosts included the volume 
work done by the expanding radio lobe. This part of the ghosts energy, which is j of its total 
energy budget, is located in the IGM in the form of heat, kinetic, and gravitational energy. 



1 During structure formation the gas follows the energy equation 

d fi n 7 ^gas \ _ d^grav 

-T, T^gas + 7 1" *grav -— , (6) 

dt y 2 7-1 e gas J dt 

if heating and cooling processes can be neglected. This equation can be integrated, if we approximate the 
potential to be time independent, yielding 

\v 2 + — Y — ^ = -(-t-grav - <t>grav,o). (7) 

Whenever the kinetic energy is small a unique relation between the potential and the specific enthalpy is 
given. 
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The normalization factor can therefore be expressed as 

gh K*r(2 — CKgh) -E'min/'-E'*) ' ^ 

where T is the incomplete Gamma-function: 

/>oo 

r(o,x) := / dtt^e- 1 . (15) 
The total number density in ghosts reads 

Finally we assume that the energy density in ghosts scales with the galaxy density to 
some power 6 g h 

£gh(n g ai) = e g h,o(n-gal/n-gal,o) bgh - (17) 

6 g h = 1 would correspond to constant energy output per galaxy. This might be unrealistic, 
since the older, and bigger galaxies in the denser regions are expected to have produced more 
radio plasma during their lifetime. We therefore also adopt 6 g h = 2 in order to have a ghost 
distribution more clustered than the galaxies. 

On the other hand, buoyant motion of the ghost might have allowed a significant fraction 
to escape from clusters and filaments of galaxies. This could have let to a distribution of 
ghosts which is less clustered than the galaxies. We mimic this case by also using b g h = 1/2. 

2.2 Guessing the Parameters 

The parameters, which are still undetermined, are e g h,o: c* g h, E m [ n , and -E*(n ga i). 

The progenitors of ghosts are RGs. We therefore try to use the observed radio luminosity 
function (RLF) of radio galaxies to get an idea about the shape and normalization of the 
ghost size distribution. The present day spatially averaged number density of ghosts is simply 
given by the time integral of their birth rate 

n-gh(-Egh) = J dth gh (E gh ,t). (18) 

The birth rate is given by the number density of visible radio galaxies divided by their typical 
lifetime ^visible- We assume that the time the galaxy is active ^active and the time the galaxy 
is visible ^visible are the same for all galaxies (i a ctive = ^visible = t Ig — 3 ■ 10 7 yr). This is a 
typical value out of t he observed range of ages of large expanded radio lobes, which are of 



the order of 100 Myr (Alexander and Leahy, 1987 ; Komissarov and Gubanov, 1994 ; Vcnturi 



et al., 1998; Slee and Roy, 1998|). Although these fairly robust age estimates are based on 



interpreting the observed radio spectral steepening in terms of the inverse Compton losses 
of the relativisitic electrons against the ubiquitous cosmic microwave background photons 
uncertainties and scatter are likely up to one order of magnitude. We note, that the assumed 
age only affects the distribution of the total ghost energy budget onto small and large ghosts 
sizes, but not the budget itself. 

We further assume that the jetpower g Je t and radio luminosity L u of a RG does not change 
during its active phase. These are very coarse assumptions. But we hope they are justified 
since we use them only to get a feeling for the ghosts size distribution. Radio galaxies with 
jetpower gj et produce therefore ghosts with the energy 

^gii = 2 q i ct tr s ( 19 ) 



with a birth rate of 



rcghl-C'gh,*) - 2 — 4 n > ( 2U ) 
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RLF b v e gh /(10 66 ergGpc- 3 ) n gh /(10 8 Gpc~ 3 ) 



PLE 0.7 2.80 0.28 

PLE 0.82 2.19 0.28 

PLE 1.0 1.91 0.28 

RLF2 0.7 5.05 1.11 

RLF2 0.82 3.32 1.11 

RLF2 1.0 2.55 1.11 



Table 1: Cosmological energy output of radio galaxies for two radio luminosity functions of 
Dunlop & Peacock Q1990Q which are used in Fig. |IJ The last column gives the total number 
density of ghosts in the ernergy range plotted in Fig. |l]. 



where n^ et {q- iet , t) dqj et is the jetpower distribution function of RGs. The number two appears 
since RG have typically two radio lobes. 

The jetpower distribution function can be derived from the radio luminosity function of 
RG with the additional assumption that there is a unique relation between radio luminosity 
and jetpower. Enfilin et al. (1997) have fit a power-law relation to the jetpower derived 
by Rawlings and Saunders (1991) of a sample of radio galaxies. Their values are based on 
minimum energy arguments and age estimates. The real energy of the radio plasma can 
easily be much higher than the minimal energy estimate by some factor / power > 1 due to the 
presence of relativistic protons, low energy electrons or deviations from equipartition between 
particle and field energy densities. A rough estimate of / power can be derived from obser- 
vations of radio lobes embedded in the intra-cluster medium (ICM) of clusters of galaxies. 
They show a discrepancy of the thermal ICM-pressure to the pressure in the radio plasma 
following minimal energy arguments by a factor of 5 — 10, even if projection effects are taken 



into account (Feretti et al., 1992). Since also a filling factor smaller than unity of the radio 
plasma in the radio lobes can mimic a higher energy density we chose /power = 3 in order to 
be conservative. 

The jetpower-radio luminosity correlation at v — 2.7 GHz then reads 

q jct = a u (L„/(Watt Hz" 1 Z^ 2 )) 6 " / powor (21) 



for which K = 0.82 ±0.07, a v = l0 45 - 28 - 26 - 22 b - ±018 ergs" 1 /i^ 2 flEnfilin et al., 1997|) . This 
relation allows us to translate the observed radio luminosity function into a jetpower distri- 
bution function. We adopt different radio-luminosity functions parameterized by Dunlop & 



Peacock ( 1990 ) and integrate Eq. [L8| using Eqs. 20 and |2lj in an Einstein de Sitter Cosmol- 
ogy up to redshift z — 10 (the result depends only weakly on the upper cutoff due to the 
strong decline of the used RLFs beyond z = 3, the pure luminosity evolution (PLE) RLF 
is only integrated up to z = 4). We further calculate the ghost distribution function for 
b u = 0.7,0.82,1 in order to show the dependence on the uncertainties. The results can be 
seen in Fig. |l|. 

Fig. |l| shows that a g h ~ 2, which we adopt in the following. A reasonable lower cutoff 
is then E ni i n — 6 ■ 10 57 erg (corresponding to L v = 10 23 Watt Hz -1 ). The ghost distribution 
does not need to have a real cutoff there, but the part which is accessible to radio observations 
of the parent radio galaxies ends there. In order to conservative we do not extrapolate to 
lower energies. 

The radio luminosity function itself becomes dominated by starburst galaxies at smaller 
radio luminosities. This gives a second contribution to the magnetization of the IGM from 
magnetized galactic winds. An extrapolation of the strongly redshift dependent, recent star- 
burst history of galaxies to high rcdshifts [z = 10) indicates that a substantial cosmological 



volume might be magnetized by such winds (Kronberg et al., 1999). However, the fields 
resulting from this process are expected to be tangled on small, galactic scales, and therefore 
less effective scattering centers for UHECR as radio ghosts. 

An integration of the ghost energy distribution functions shown in Fig. [l| gives the total 
energy output of radio galaxies per volume. The results do not depend on the assumed 
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lifetime £ rg of radio galaxies, they scale linearly with / power , and depend only moderately on 
b v . Detailed results are given in Tab. |[ The energy input in form of radio plasma is roughly 
e g h = 10 66 -" 67 ergGpc~ 3 , it is distributed over n g h = 10 7 - 8 ghosts per Gpc 3 , which have an 
average energy of e g h = 10 58 — 59 erg. Note that the X-ray background, which is believed to 
be dominated by AGN emission, corresponds to an injection energy of w 3 ■ 10 67 erg Gpc~ 3 
(Chokshi & Turner 1992; Soltan 1982). Either the X-ray energy losses of AGNs exceed the 
radio plasma release, or the jetpower is strongly underestimated here and / powcr = 10. ..30 
would be more realistic. Even such a strong energy input into the IGM seems to be consistent 



with present the day limit to the Comptonizarion of the IGM (Enfilin and Kaiser, 200C ; EnBlin 
et al., 1998b| ) 



The next parameter we have to fix for our model is the maximal ghost energy. E*(n sa .i) 
can be roughly estimated from observations as shown in the following. The progenitor of large 
ghosts in rare regions should be giant radio galaxies (GRGs). Their linear size distribution 
seems to have a sharp cutoff at 3 Mpc (Ishwara-Chandra & Saikia 1999). There is only 
one GRG significantly exceeding this: 3C236 has a linear size of 5.7 Mpc h^. This could 
be - for example - due to a very low pressure in the environment of 3C236. For a typical 
galaxy filament environment, we assume the linear size of GRG to have a cutoff at 3 Mpc 
h^Q . The diameter of one of the two radio lobes is roughly -| of the linear size of the whole 
GRG. We adopt therefore a maximal radius of our (spherically approximated) ghosts in the 
dilute environment of galaxy filaments of r* = 0.5 Mpc h§Q . The distribution function of 
ghost sizes 

3(r/r») 2 " 3 "^ / r 3 \dr 

fo(r) dr = — '- — — exp — ^ — (22) 

has indeed a very sharp cutoff at r*, as required by the observation of GRG. This r* is related 
to E„ via Eqs. || and g: 

E* = f nr^P^f *rle B = ^rlB% (23) 

Equipartition (or minimum energy) magnetic field strength of GRG are of the order of a few 
/iG (Ishwara-Chandra & Saikia 1999), so that we adopt a conservative value of £? g h = 1/iG 
for ghosts located in galaxy filaments. From this we find that E* = 1.6 • 10 60 erg, and that 
a corresponding RG had ejected twice as much energy (3.2 • 10 60 erg), which seems to be 
moderate. Comparison with Fig. [I] shows that ghosts with 10 — 100 times higher energy are 
expected from the radio counts (if t rg is as big as 3 • 10 7 yr). Such ghosts are also required 
since a 10 60 erg ghosts in a dense cluster environment (n e = 10~ 3 cm~ 3 — > £> g h = 25/iG) has 
a diameter of 0.1 Mpc, which is clearly below the observed maximal size of radio lobes inside 
clusters. In order to allow a diameter of 0.5 Mpc inside a cluster one needs E* w 10 62 erg. 
We therefore assume that the maximal energy of a ghost scales linearly with the gas density: 

* ~ 10M "MF^- 10M ^(S) 1 """' (24) 

where n ga \ iQ is the galaxy density in the center of a typical cluster of galaxies. Such a scaling 
might be justified due to the fact that the most massive galaxies, which probably host the 
most powerful AGNs, are the central elliptical galaxies within clusters. 

For the purpose of our simulation it also makes sense to define the average total cross 
section for the assemble of ghosts: 

/j CI \ 2 2 -T(| — Ckgh, £'miii/£'*) 

drf{rUr =nr * F(l-a gh ,iW^) ^ 

2.3 Constructing the Ghost Distribution 

The scaling relations in Eqs. |l0|, [IT], and |24| are normalized for a reference density equal to 
the central galaxy density of a galaxy cluster. We estimate n ga i j0 by identifying the peaks 
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in our galaxy distribution. The energy density of ghosts at our reference gas density n ga i j0 
(inside clusters) is given by 



e K h,o 



e gh V 



J dV (n gai /n gai:0 ) 1 



(26) 



where e g h = 10 66 " 67 erg Gpc~ 3 and V is the volume of our galaxy catalog. Whenever we use 
the CfA redshift catalog we exclude volumes close to the galactic plane in our estimate of V 
and the above integral due to the under-sampling of galaxies there. 

If the ghosts would get disrupted and mix with the ambient gaseous medium a signifi- 
cant magnetization of the IGM would result. Assuming complete mixing, we estimated the 
resulting magnetic field strengh as a function of the local cosmic matter density in a few 
scenarios (e g h = 10 66 ' " 67 ergGpc~ 3 , 6 g h = 0.5 and b g h = 1). The field strength as a function 
of density is plotted in Fig. and the volume fraction of the different density regimes in 
Fig. Such magnetic fields are not excluded by Faraday rotation measurements, which 
give limits of -Bigmf < nG for a homogeneous background field with Mpc coherence-scales 
( Kronberg, 1994b ), or -Bigmf < ^G if cosmologi cal fields are concentrated in the larg e scale 
structures like sheets and filaments of galaxies ( Ryu et al., 1998a ; Blasi et al., 1999] ). The 
implied values for cluster of galaxies of -Bicm ~ 0.3...3/xG are also consistent with Faraday 
rotation measurements of magnetic fields. Note, that shear flows due to cosmic structure 
formation can signific antly increase the strength of IGM magnetic fields ( Dolag et al., 199S ; 
Rocttigcr et al., 1999|). 



2.4 The Radio Ghosts Produced by Our Galaxy 

The formation of the Mbh = 2.6 • 10 6 M & black hole in the center of our galaxy (Eckart et 
al. 1997) was probably accompanied by dissipative processes as radiation and ejection of 
relativistic plasma in a temporary AGN. We assume that a fraction ediss = 0.1 of the rest 
mass energy of the accreted matter was dissipated and that e p i aS ma = 0.5 of this energy went 
into the relativistic plasma. The radio ghosts formed by the remnant of this radio plasma 
got therefore an energy budget of 

E gh = epla ^ ma£diss M bh c 2 = 2.6 • 10 59 erg . (27) 

1 ^diss 

The term 1 — ediss in the denominator corrects for the mass discrepancy between accreted 
mass and final black hole mass due to dissipation of energy. Eqs. [j], ||, and |^ imply that the 
filled volume is given by 

^ = ^ = a.083Mpc*G§)-\ < 28 > 

Comparing this to the gyroradius r g of a 10 20 eV cosmic ray proton shows the possible 
importance of the ghosts of our own galaxy: V^h ~ 66 r 3 (-B//iG). 

But it cannot be assumed that the ghosts produced by our own galaxy are located in the 
direct vicinity of our galaxy since the observed peculiar velocity of 162 km/s with respect to 
our own galaxy group (Rauzy & Gurzadyan 1998) indicates differential motion between the 
Milky Way and the local IGM. Radio ghosts produced by other members of the local group 
may be nearby today, but this is already included into our statistical approach described in 
Sec. ||, where the density of ghosts scales with the number density of galaxies. 



3 UHECR Flux at Earth 



3.1 Scattering by Radio Ghosts 

Equation (5) can be used to calculate the intensity of the magnetic field inside ghosts. The 
estimation of the topology of the internal field is, however, a much more difficult matter. 
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The generation of any ghost was a turbulent event fed by the parent radio galaxy. This 
process certainly produced vortices at all scales. Ghosts, however, have been dynamically 
decoupled afterwards and, therefore, smaller scale wiggles should have been dissipated over 
time. Hence, most of the remaining power in the field probably lays at large scales at present. 

In Fig. [2] we exemplify one possible internal magnetic field topology. The shown magnetic 
field has null perpendicular component at the surface and is therefore decoupled from the 
environment field. The amplitude of the fluctuations obeys a power law power spectrum of 
spectral index £ = 4. 

In order to analyze the scattering capability of ghosts and to estimate how critical is the 
uncertainty in the actual topology of the field to the present study, we performed numerical 
simulations using spherical radio ghosts, similar to the one in Fig. |. Different power spectral 
indexes £ and radius were considered. We illuminate the ghosts with a beam of test particles, 
protons, with a power law energy spectrum, dN/dE cx E~ 3 , and energy E > 4 x 10 19 eV. The 
angular distribution of the scattered particles is shown in Fig. 3. a for ghosts of different sizes 
and £ = 5/3, 3 and 4. The radii of the ghosts are given in terms of the average gyroradius of 
the injected UHECR spectrum. As a comparison, an isotropic scattering particle distribution 
is also shown in the same figure. It can be seen that ghosts are very efficient at scattering 
UHECR regardless of power spectral index. Forward scattering becomes dominant only at 
the low end of the ghost size distribution, when the ghost's radius becomes comparable 
to the average gyroradius of the particles. Figure 3.b is the same as figure 3. a but for a 
monoenergetic proton injection at E = 10 20 eV. The same ghosts radii as in fig. 3. a (physical 
units) were used, and so they are down by a factor of two in units of UHECR gyroradius 
in fig. 3.b. The stronger dominance of forward scattering for the smallest ghosts can be 
appreciated, as well as for the flattest fluctuation spectrum (£ = 5/3). Nevertheless, as 
previously argued, small scale magnetic turbulence inside a typical ghost should have been 
dissipated over time, and a steep spectrum of magnetic fluctuations is expected. 

Therefore, given the current uncertainties in radio ghost shape, size, magnetic field topol- 
ogy and spatial distribution, isotropic scattering is, very likely, an acceptable assumption and 
it will be used in what follows. 



3.2 Numerical Model 

We use Monte Carlo numerical simulations to track UHECR propagation through the inter- 
galactic medium and to evaluate their arrival distribution at Earth. 

We start from the basic assumption that most UHECR are protons of extragalactic origin 
whose sources aggregate spatially as either luminous or cold DM. 

Charged particles, even at the extreme energies considered, are coupled to the intervening 
magnetic fields. From this point of view, the propagation region can be divided in several 
components: (a) the sources and their immediate neighborhood, (b) the intergalactic medium, 
(c) the galactic halo, (d) the galactic disk, (e) the heliosphere and (f) the magnetosphere. 

The close environment of the sources can be included inside the sources just by re-defining 
them as large as needed. This should be acceptable in general since, the newly defined source, 
would still be below the resolution capability of any current experiment. 

The heliospheric and magnetospheric fields, despite their relatively large values, arc con- 
fined to small volumes and, from the point of view of deflection, are negligible. 

The galactic disk magnetic field is one to two orders of magnitude smaller than the latter 
fields, but fills regions at least 8 orders of magnitude larger. However, exception made of 
lines of sight crossing directly the galactic bulge or a relatively narrow strip surrounding the 
galactic plane, the disk field has probably little impact on UHECR deflection. 

Consequently, the two main factors responsible for blurring pointing information, are the 
intergalactic (IGMF) and the galactic halo magnetic field. 

Few is known regarding the large scale structure of the IGMF and only sparce rotation 



measurements are available to constrain its intensity ( Kronberg, 1994b ; Vallee, 1997 ). Basi- 
cally, two extreme scenarios can be imagined that satisfy the observational constraints. In 
one of them, the field is compressed inside walls and filaments, where it attains high intensity 
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(~ 0.1 — 1/iG) and high de gree of correlation (over scales of up to tens of Mpc), leaving 
unmagnetized voids behind ( Ryu ct al., 19981: ). We will call this the laminar IGMF model. 
A second possibility is that of a more evenly distributed field, whose intensity scales some- 
how with matter and is coherent on scales L c w 1 Mpc (B /nG)~ 2 . A numerical model for 
the latter case is an ensemble of cells of characteristic size L c , with uniform field inside but 
randomly oriented with respect to adjacent cells. The galaxy density is used to scale B and 
L c . We call this model the cellular IGMF. 

The first IGMF scenario can affect propagation in a radical way, with UHECR mainly 
constrained to the interior of large scale structures and drifting along walls and filaments. 
The observable flux depends critically on the location of the observer and on the details of 
the actual topology of the field inside the GZK-sphere. Isotropy can be attained in a rather 
natural way, but the results are strongly model dependent. The propagation of UHECR 
in this IGMF has been treated in an earlier work (Medina Tanco, 19981), and will not be 
considered here. 

A similar uncertainty exists around the topology of the magnetic field inside the galactic 
halo. Nevertheless, unless there is a large scale magnetized galactic wind that structures 
this field and maintains a large azimuthal component up to large galactocentric distances, 
the halo should be unimportant regarding isotropization. If a magnetized galactic wind 
does exist on the scale of several hundreds of kpc , then it would b e the dominant factor i n 
the determination of the UHECR flux at Earth ( Ahn ct al., 200C ; Bicrmann et al., 200C ); 
therefore, extragalactic propagation details would be mostly irrelevant. 

Therefore, we restrict the present analysis to a cellular IGMF and neglect the galactic 

halo. 

The same procedure as in ( Medina Tanco, 1997 ; Medina Tanco, 1998a ) is used in the 
description of the cell-like spatial structure of the IGMF. The cell size is given by the corre- 
lation length, L c cx Bj2 MF {r). The intensity of the IGMF, in turn, scales with gas density 



as Bjgmf oc n^ as (r) and the proposed IGMF value at the Virgo cluster (~ 10 7 G, (Arp 



1988)) is used as the normalization condition. We use rj = 0.5, as a compromise between a 



frozen-in field (r) — 2/3) and Vallee's (1997) estimate (r) ~ 0.35), which may be too flat due 
to the assumed values for the magnetic field in superclusters and larger scales. Nevertheless, 
tests have been conducted for different values of rj covering the previous interval, and the 
scaling is not critical to our conclusions. 

The formulation of section 2 is used to define the size and spatial distribution of radio 
ghosts. This requires the knowledge of the gas density distribution inside the simulation 
volume. Actually, as galaxies are easier to survey over large volumes, equation(9) is used 
to transform between galaxy and gas density distributions. There are, however, serious an 
unavoidable bias and sampling problems inherent to galaxy surveys. Therefore, we relay on 
both galaxy surveys and cold DM large scale structure simulations to perform independent 
evaluations. 

The sources of UHECR are distributed according to either the galaxy or cold DM distri- 
butions respectively. 

The 1999 version of the CfA catalog ( Huchra et al., 1992 ) is used to characterize the 
galaxy distribution. 

Cold DM simulation data are from Springel et al. (in preparation). They carried out 
simulations that mimic the Local Universe. The initial conditions of these simulations have 
been constrained by the redshift survey of IRAS galaxies. As a result, the simulations develop 
the same local large-scale structure (e.g., the Great Attractor and Cetus Wall; clusters like 
Virgo and Coma are also found at the right place). 

The scenario is completed by the introduction of radio ghosts according to the formulation 
of section 2. 

Relativistic test particles (UHECR protons) are injected at the sources with a spectrum 
dN/dE cx E~ 2 and propagated through the intergalactic magnetic field up to the detector 
on Earth. Adiabatic energy losses due to redshift, pair production and photo-pion produc- 
tion due to interactions with the cosmic microwave background radiation (CMBR) are also 
included (Bcrczinsky and Grigor'cva, 198? ; Achterberg et al., 1999| ). 
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The flux at Earth can be divided into two components: (a) a direct radiation field, 
constituted by particles that fly from source to detector without encountering ghosts and (b) 
a diffuse radiation field, comprising particles which underwent at least one encounter with 
a radio ghost. 



3.3 Different Scenarios 

Fig. |] is meant as a control. It shows the Aitoff projection of the two-dimensional arrival 
probability density (galactic coordinates with the anti-galactic center at the center of the 
figure), for s ources distributed according to nearby luminous matter (CfA catalog, ([Huchra 



ct al., 1992 )) inside 100 Mpc and no radio ghosts. The same procedure as in ( Medina 



Tanco, 1997 ; [Medina Tanco, 1998a ) is used in the description of the intergalactic magnetic 



field (IGMF): a cell-like spatial structure, with cell size given by the correlation length, 
L c oc BfQ MF {r). The intensity of the IGMF, in turn, scales with luminous matter density, 



Pgai as Big mf oc p®gi{ r) flVallee, 1997[ ) and the observed IGMF value at the Virgo cluster 
(~ 10~ 7 G, ( Arp, 1988 )) is used as the normalization condition. The mask covers the plane 
of the galaxy, where the actual distribution of galaxies is not well known due to obscuration 
by dust. The curved, thick line is the celestial equator. Northern hemisphere is the sky patch 
to the right, enclosed by that line. 

Superimposed on the figure are the available events with E > 4 x 10 19 eV observed by 



AG ASA (47 events (|Takeda, 1999| )), Haverah Park (27 (|Reid and Watson, 1980| )), Yakutsk 



(24 (Afanasiev, 1995)) and Volcano Ranch(6 (Linsley, 1980)) 



The arrival probability contours trace roughly the local large scale structure. Distinguish- 
able observational signatures should be expected towards the region of the Southern branch 
of the supergalactic plane (to be observed in the near future by the Auger experiment) at 
/ ~ 45°, the lines of sight to the more distant Pisces-Perseus wall and Perseus cluster and, 
very prominently, towards a large area surrounding the Virgo and Ursa Major clusters. 

The actually observed distribution of UHECR is clearly much more isotropic than what 
one would expect under the implicit assumptions in Fig. 0. 

In figures 5, 6 and 7, we show how the UHECR arrival distribution function would be 
modified by the presence of radio ghosts in the intergalactic medium for different values of 
b g h and a g h = 2.0. The results are insensitive to a g u in the range of interest. The most 
important parameter is bgh, defined in equation (17), which tells how are ghosts distributed 
with respect to galaxies. b g h = 1 implies the same spatial distribution as galaxies and larger 
values a more clustered distribution. 

It can be seen from Fig. || and [| that clustered distributions of ghosts produce no 
noticeable effects in the observed UHECR flux at Earth. 

In fact, if these angular distributions were observed with the same exposure in declination 
as AGASA's aft er 7.5 years of integration, the observed (ensemble averaged) amplitude of the 
first harmonic ( Linsley, 1975] ) would vary only by ~ 5% when going from Fig. | to |. This 



is basically because a distribution of scatterers that is more clustered than the sources is, 
in practical terms, equivalent to increase the physical size of the sources, which is irrelevant 
from the point of view of all sky anisotropies. 

The result is different when the scatterers are distributed in a larger volume than the 
sources. This is the case in figure 7, where b g h = 0.5. In this scenario it is assumed that 
e.g. radio ghosts buoy out of cosmological structures, creating thick halos around walls 
and filaments, permeating voids to some extent. This diminishes considerably the direct 
component (Fig. 0a) and accounts for a considerable increase in the diffuse component (Fig. 
[7(3) which becomes, by far, dominant. The composite flux (Fig. ^c) still shows a smooth, 
large scale gradient towards the region of the Virgo cluster but is much more isotropic that 
any of the previous scenarios. 

A more quantitative picture can be obtained from Fig. ||. An ensemble of 1000 indepen- 
dent samples, of 50 UHECR each, was built by observing the arrival distribution function in 
Fig. 0c with an exposure equivalent to AGASA's. The resultant distribution function of the 
amplitude of the first harmonic is shown, in comparison with the actually observed values 
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by Y akutsk, AGASA, Haverah Park and Volcano Ranch (see, ( Medina Tanco and Watson 



1999)). It can be seen that b g h — 0.5 is the upper limit for the clustering of ghosts that is 



compatible with the observed anisotropy. However, it is also clear that that an even smoother 
distribution of ghosts in intergalactic space would be much more acceptable. Increasing the 
volume average ghost energy density by a tolerable factor of 3 would help to move the results 
in the right direction, but it wouldn't be a determinant factor by itself. Rephrased in other 
words, and looking at Fig. |[ the mean-free path of UHECR for interactions with ghosts in 
the intergalactic medium, should be reduced to less than a few Mpc at any location in the 
propagation region, i.e., walls and filaments as well as voids, to obtain a level of isotropy 
consistent with observation. 

Whether such a smooth spatial distribution of radio ghosts is physically reasonable is a 
very complicated matter which we leave open. The possibility that radio ghosts are able 
to escape by buoyancy from the denser environments in which they were produced exists. 
However, specific large scale structure hydrodynamic simulations should be performed to 
decide whether buoyancy is actually operational and, furthermore, a deeper knowledge of the 
composition and mechanical properties of the radio plasma is necessary in order to decide if 
ghosts are stable enough to survive as an entity under these conditions. 

As was mentioned previously, there are uncertainties associated with unavoidable biases 
and sampling incompleteness associated with galaxy surveys, as is the case with the CfA 
catalog used up to here. This problem can be specially critical in this analysis, since the 
results depend on an absolute normalization of the density as a function of depth into the 
local universe. To check the extension of the distortions occurring in our previous analysis, 
we repeated our calculations using the distribution of cold DM, calculated by large scale 
structure hydrodynamic simulations (Springel, in preparation). They carried out simulations 
that mimic the Local Universe, developing the observed local large-scale structure. 

Fig. |l^ and [ll] show the arrival flux for a g h = 2.0 and b g h = 2.0 and 0.5 respectively 
for the simulation using the CDM model. The resultant UHECR flux is plotted in an Aitoff 
projection in supergalactic coordinates. The large spot near the center of Fig. [To] is the 
cosmic ray image from the Virgo cluster. The central band is the supergalactic plane. Note 
that the supergalactic plane is not such a neatly defined structure in UHECR with average 
energy E > 4 x 10 19 eV. It is clear that, regarding isotropization of the UHECR flux, the 
same kind of effect is present in these new scenarios: the dominant signature from Virgo at 
bgh > 1 is strongly diluted at b g h < 0.5. 



4 Discussion 



The, to some extent unexpected, isotropy in the UHECR flux observed particularly by 
Yakutsk and AGASA, is a very strong constraint to production and propagation models 
alike. Most of the scenarios proposed in the literature should produce a sizable anisotropy 
as either extragalactic luminous or dark matter is normally associated with the invoqucd 
particle sources. If that is really the case, then it is our understanding of the topology and 
intensities of the intervening magnetic fields that is critically incomplete. The problem could 
reside inside the network of walls and filaments, in the interior of the large surrounding voids, 
or even in our nearby environment, namely the Galaxy halo. 

In the present work we explore the possibility that radio ghosts, blobs of magnetized 
radio plasma remnant from past periods of activity in radio galaxies, being able to scatter 
UHECR in the intergalactic medium. Such a process could, in principle, degrade the direct 
incoming flux from the sources and build up a diffuse UHECR component large enough to 
be responsible for the observed degree of isotropy. 

Our results show that, over the most conservative region of the radio ghost parameter 
space, such isotropization is not possible. This stands not from an inability of ghosts to 
scatter UHECR, but mainly from the fact that, under general conditions, ghosts should tend 
to cluster more strongly than the sources of the particles. 

If, however, radio ghosts are able to buoy out into the surroundings of the dense large scale 
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structures and into voids, while surviving the process, UHECR isotropy could be obtained in 
those cases in which the mean free path for interactions with ghosts is reduced below some 
few Mpc all over the propagation region (<~ b g h < 0.5). 

Finally, one can question the basic assumption of this work, that radio ghosts are able 
to survive for cosmological times as intact objects. If radio ghosts mix with the ambient 
medium a significant magnetization of the IGM should result, which can be sufficiently 
strong in order to explain observed magnetic fields in clusters of galaxies. These magnetic 
fields would have short coherence length, and therefore do not violate Faraday rotation 
measurements. The ability of the fields to deflect UHECR particles would be increased 
compared to the case of intact, non- mixing ghosts for the following reasons: the average 

2 /3 

deflection of particles in a ghosts vicinity is proportional to Vgh -Bgh (cross section U gh ' 

1/3 

times scattering angle, the latter is proportional to -B g h times V^ h ). For a fixed magnetic 
energy budget E Bgh = Vgh B^ h / (8tt) of the ghosts, this implies that the scattering efficiency 
scales with V^ 2 and therefore increases with increasing size. For fixed magnetic flux 

(i?gh oc V~ 2 / 3 ) the scattering efficiency inreases with the size as V^[ Z . This demonstrates 
that relaxing the hypothesis of compact ghosts makes this scenario even more powerful. 
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Figure Captions: 



Fig. [l]. Energy distribution of radio ghosts multiplied by E^ h for better display. The lines 
are labeled with b v , the exponent of the radio luminosity-jetpower correlation. The solid 
lines are derived from the pure luminosity evolution (PLE) radio luminosity function of 
Dunlop & Peacock (1990). The dashed lines from their radio luminosity function No. 2 
from the MEAN-z data (RLF2). This latter luminosity function has an exponential cutoff at 
luminosities of 10 28 Watt Hz -1 . The curves are plotted for the range of radio luminosities 
between 10 23 to 10 29 Watt Hz -1 . A horizontal line in this diagram corresponds to a g h = 2. 
E.g. the PLE with b v = 1.0 has a low energy index of a g h = 1.75. 



Fig. g. Example of a possible internal radio ghost magnetic field topology. The vectorial 
representation of the component of the magnetic field parallel to the plane of a slice, cut 
through the middle of the ghost, is shown. The length of the arrows scales linearly with 
magnetic field intensity. The magnetic field is tangent to the surface of the ghost. The 
amplitude of the fluctuations obeys a power law power spectrum of spectral index £ = 4. 
See text for further details. 



Fig. ||. (a) Angular distribution of scattered UHECR for ghosts of different radius (in units 
of proton gyroradius) and values of £. The incident particles have a power law spectrum, 
dN/dE oc E~ 3 for E > 4 x 10 19 eV. For comparison, the thick line shows an isotropic angular 
distribution. Forward scattering becomes dominant only as the ghost radius becomes 
comparable to the UHECR gyroradius. (b) Same as (a) but for protons at a fixed energy, 
E= 10 20 eV. The same ghosts are used as in (a) and so their radii are down by a factor of 
approximately two in units of particle gyroradius. 



Fig. |j. UHECR propagation without radio ghosts. Aitoff projection of the arrival probabil- 
ity density. Galactic coordinates are used with the antigalactic center at the center of the 
figure. UHECR sources are distributed according to nearby luminous matter (CfA catalog). 
Data points correspond to cosmic rays observed by AGASA, Haverah Park, Yakutsk and 
Volcano Ranch with E > 4 x 10 19 eV. Clearly, the observations are more isotropic than 
what should be expected from the model. 



Fig. H UHECR received at Earth for a g h — 2.0, b g h = 2.0 and sources distributed as 
luminous matter. 



Fig. H UHECR received at Earth for a g h = 2.0, b g h = 1.0 and sources distributed as 
luminous matter. 



Fig. [7] UHECR received at Earth for a g h = 2.0, b g h — 0.5 and sources distributed as 
luminous matter, (a) UHECR that encounter no ghost from source to Earth; (b) UHECR 
that underwent at least one scattering off a RG. Comparing this figure with Fig. |^, || and 
pi it can be seen that ghosts must have a more extended distribution than the sources 
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(bgh < 0.5) in order to isotropize the UHECR flux at Earth. A tolerable (3-fold) increase in 
average ghosts energy density goes in the right direction but is not a determinant factor in 
itself. 



Fig. | Amplitude of the first harmonic for UHECR with E>4x 10 19 eV for a = 2.0 and 
b g h = 0.5. 1000 independent samples of 50 UHECR each were selected from the distribution 
in Fig. |^c with the same sensitivity in declination as AGASA's. The amplitudes of the 
first harmonic obtained for these samples arc shown in the histogram. As a comparison, 
the amplitudes obtained by AGASA, Haverah Park, Yakutsk, and Volcano Ranch (Medina 
Tanco and Watson, 199E) are also shown. It can be seen that, in order to isotropize the 



UHECR flux, the condition b g h<0.5 must be fulfilled. That is, ghosts should be much more 
smoothly distributed in space than the galaxies that originated them. 



Fig. ^ UHECR mean free path [Mpc] as a function of 6 g h, the bias between galaxies and 
ghost distributions, and gas density relative to the reference level. Ghosts are able to 
isotropize the UHECR flux only for b g h < 0.5, for which the mean free path is smaller than 
a few Mpc both at voids and structures. 



Fig. UHECR received at Earth for a g h — 2.0, b g h = 2.0 and sources distributed as 
ACDM. Supergalactic coordinates are used in this figure. The strong excess in the arrival 
direction distribution corresponds to the Virgo cluster. 



Fig. |TT] UHECR received at Earth for a g h — 2.0, b g h — 0.5 and sources distributed as 
ACDM. Again, isotropization is obtained only for a ghost population that is much more 
extended than the UHECR sources. Supergalactic coordinates are used in this figure as in 
Fig. 0. 



Fig. [12] Distribution of injected magnetic field strength as a function of background density, 
in the case that all radio ghosts completely mix with the ambient medium. The volume 
fractions of the different density regimes can be seen in Fig. |l3|. 



Fig. The volume fractions of the different density regimes (dV/d\og(n gas /n gas fi)/ J dV) 
computed from the simulated dark matter distribution. The median density is marked. A 
count in cell box with size of 1.5 Mpc was used. 
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